REPORT  DOCUMENTATION  PAGE 


♦ 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0180),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 


FINAL 


4.  TITLE  AND  SUBTITLE 

FROM  DYNAMICAL  SYSTEMS  TO  COHERENT  STRUCTURES  SHADOWING 
CHAOS  WITHIN  TURBULENCE 


01  MAR  93  TO  29  FEB  96 


5.  FUNDING  NUMBERS 

F49620-93-1-0172 

61102F  2304/CS  AND  2304/AS 


6.  AUTHOR(S) 

DR  BASEL  NICOLAENKO 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
ARIZONA  STATE  UNIVERSITY 
BOX  871603 

TEMPE,  ARIZONA  85287-8013 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFOSR/NM 

1 1 0  DUNCAN  AVE  SUITE  B 1 1 5 
BOLLING  AFB  DC  20332-8080 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

<H  U-OL-iSO 


1 10.  SPONSORING/MONITORING 
I  AGENCY  REPORT  NUMBER 

F49620-93-1-0172 


11.  SUPPLEMENTARY  NOTES 


1 9960822  233 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 
APPROVED  FOR  PUBLIC  RELEASE: 
DISTRIBUTION  UNLIMITED 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

SEE  REPORT 


mio  qUAUTS  fflSFBttriB)  « 


114.  SUBJECT  TERMS 


115.  NUMBER  OF  PAGES 


116.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  Form  298iRev.  2-89)  (EG) 

Prescribed  by  ANSI  ^d.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


Shadowing  Chaos  Within  Turbulence 
Grant  AFOSR  F49620-93-10172 
Final  Technical  Report 

Principal  Investigators: 

Basil  Nicolaenko  and  Alex  Mahalov 
Department  of  Mathematics 
Arizona  State  University 
Tempe,  AZ  85287-1804 

e-mail:  byn@stokes.la.asu.edu,  alex@taylor.la.asu.edu 

Abstract 

The  first  focus  is  on  parametrization  and  decomposition  of  turbulence  in  geophysical 
flows.  Such  flows  are  described  by  3D  Euler/Navier-Stokes  equations  with  Boussinesq 
stratification  and  uniform  rotation.  Within  the  context  of  rotation  in  thin  atmospheric 
layers,  we  demonstrate  that  geophysical  turbulence  splits  between  2D  turbulence  for  the 
geostrophic  component  and  phase  turbulence  for  the  ageostrophic  component.  By  solv¬ 
ing  low-dimensional  dynamical  systems,  we  obtain  exact  formulas  for  phase  turbulence 
and  phase  scrambling  in  atmospheric  dynamics.  Numerically,  this  leads  to  rigorous  non¬ 
stiff  operator  splitting  algorithms  for  geophysical  flows.  The  second  focus  is  on  inertial 
stability  of  such  algorithms.  The  new  concept  of  inertial  stability  guarantees  that  long 
time  simulation  of  turbulence  via  scientific  computing  is  indeed  statistically  relevant  to 
real  physical  turbulence. 
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1  Research  Accomplishments 


1.1  Global  Regularity  and  Integrability  of  3D  Euler  and  Navier-Stokes 
Equations  for  Uniformly  Rotating  Fluids 

There  are  many  important  engineering  and  geophysical  problems  in  which  rotation  signifi¬ 
cantly  modifies  properties  of  fluid  flows.  In  particular,  large-scale  atmospheric  and  oceanic 
flows  are  dominated  by  rotational  effects  and  the  impact  of  shallowness  on  the  large  scale 
dynamics.  These  problems  attracted  attentions  of  many  physicists  and  mathematicians 
beginning  with  works  of  Laplace  [1775-1779],  Kelvin  [1880],  Poincare  [1910]  (see  Greenspan 
[1968]  and  Pedlosky  [1987]  for  more  references).  Long-time  computation  of  geophysical 
flows  using  unmodified  Euler  equations  is  prohibitive  due  to  severe  accuracy  and  time  step 
restrictions.  The  disparity  of  time  scales  leads  to  problems  in  the  numerical  solution  of 
the  equations  because  the  Courant  number  is  determined  by  the  fastest  time  scale  l/fio 
(fio-  angular  velocity  of  the  background  rotation)  and  therefore  limits  the  time  step  which 
makes  explicit  solution  impractical  (Browning,  Holland,  Kreiss  and  Worley  [1990]).  One  of 
the  major  difficulties  encountered  is  the  influence  of  the  oscillations  generated  by  the  rota¬ 
tion  (see,  for  instance  Bourgeois  and  Beale  [1994],  Lions,  Temam  and  Wang  [1994],  Babin, 
Mahalov  and  Nicolaenko  [1995a]  and  [1995b],  Chemin  [1995],  Grenier  [1995a]  and  [1995b] 
and  more  references  in  the  cited  papers,  see  also  the  list  of  references).  We  develop  here 
the  techniques  which  allow  us  to  overcome  this  difficulty  and  in  fact,  the  oscillations  help 
to  solve  the  problem  explicitly. 

We  consider  flows  uniformly  rotating  with  a  constant  angular  velocity  SIq  about  the 
13  axis.  Let  Lq  be  the  characteristic  horizontal  length  scale  and  Uo  be  the  characteristic 
velocity.  The  problem  is  made  dimensionless  by  referring  lengths  to  Lq,  velocities  to  Uo  and 
time  to  LoIUq,  which  leads  to  the  introduction  of  a  dimensionless  rotation  rate  SI  {Ro  —  1/Sl 
is  the  Rossby  number),  defined  by  fl  = 

In  a  frame  of  reference  rotating  with  constant  angular  velocity  Si  about  the  13  axis,  the 
inviscid  3D  Euler  equations  in  the  dimensionless  variables  have  the  form 


StU  4-  U  •  VU  +  2flJU 


-Vp, 


(1) 


v-u 


0. 


Here  U  =  is  the  velocity  field,  J  is  the  rotation  matrix,  2fiJU=2fie3  x  U  is 

the  Coriolis  force  term  and  p  is  a  modified  pressure.  We  consider  solutions  of  Eq.  (1)  27raj- 
periodic  in  every  Xj  taken  along  Cartesian  axes  j  =  1,2,3;  ai  =  1,0  <  03  <  a2ij0  < 
^3  <  ^31-  The  linear  version  of  (1)  was  extensively  studied  by  Poincare  in  1910  with  the 
introduction  of  Poincare’s  velocity  variables  and  Poincare’s  inertial  waves.  Poincare’s  waves 
are  well-known  in  the  geophysical  literature  (e.g.  Greenspan  [1968]). 

In  our  previous  work  (Babin,  Mahalov  and  Nicolaenko  [1995a])  we  have  used  Poincare 
velocity  variables  and  averaging  Euler  equations  in  time,  derived  new  extended  2D/reduced 
3D  Euler/Navier-Stokes  equations  for  the  Poincare’s  velocity  variables.  The  structure  of 
these  equations  is  determined  by  the  structure  of  resonant  cones  in  the  Fourier  space  of 
wavenumbers;  after  time  averaging  over  the  fast  time  scale  1/0,  these  resonant  cones  gen¬ 
erate  a  non-trivial  limit  to  the  Euler  equations  as  0  — >  +00.  The  method  is  closely  related 
to  time-averaging  methods  in  celestial  mechanics  (Poincare  [1892],  Moser  [1973],  Siegel  & 
Moser  [1971]).  In  this  work  we  resolve  the  small  divisor  problem  for  resonances  in  rotating 
3D  Euler  equations.  The  resonant  set  in  the  Fourier  space  of  wavenumbers  is  the  inter¬ 
section  of  five-dimensional  cones  with  a  lattice  in  which  depends  on  a2,a3;  we  call  this 
set  resonant  cones.  Through  a  remarkable  factorization  of  the  corresponding  small  divi¬ 
sor,  we  prove  that  the  resonant  cones  in  Fourier  space  further  reduce  to  the  union  of  four 
three-dimensional  cones  if  the  spatial  domain  itself  is  non-resonant;  that  is  the  ratios  02,^13 
of  lengths  of  edges  of  the  period  parallelepiped  are  non-resonant  {corresponding  resonant 
domain  ratios  l/a^jl/a^  are  of  zero  Lebesgue  measure  in  R,  we  assume  everywhere  below 
that  such  domain  ratios  are  non-resonant).  Therefore  aspect  ratios  are  non-resonant  with 
probability  one.  An  important  case  of  non-resonant  domains  is  small  03  (thin  domains). 
We  derive  uniform  estimates  for  the  small  divisor  in  the  complement  of  resonant  cones,  for 
non-resonant  domains.  These  are  essential  for  the  sharp  estimates  of  the  error  (in  terms  of 
the  anisotropic  Rossby  number)  in  our  strong  convergence  theorems.  Without  such  small 
divisor  estimates,  one  can  only  state  non-uniform  convergence  theorems  with  an  error  esti¬ 
mate  0(1),  such  as  done  in  Grenier  [1995a]  and  [1995b].  The  resolution  of  the  small  divisor 
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problem  is  at  the  very  heart  of  our  global  existence  theorems  for  3D  rotating  Euler  and 
Navier-Stokes  equations  presented  in  this  work. 

In  this  work  we  study  small  anisotropic  Rossby  number  flows.  The  anisotropic  Rossby 
number  Roa  is  defined  as; 

Roa  =  azIQ.  =  a^Ro.  (2) 

The  vertical  aspect  ratio  in  general  is  not  assumed  to  be  small;  as  a  matter  of  fact,  its 
smallness  improves  error  estimates  in  the  convergence  theorems.  The  anisotropic  number 
is  small  in  the  following  four  important  situations: 

(a)  rapidly  rotating  flows  (fi  >>  1,  but  aj  =  0(1)). 

(1)  rotating  flows  in  shallow  layers  (03  <<  1  but  Q  =  0(1)). 

(c)  intermediate;  Roa  =  03/0  =  a^Ro  <<  1. 

(d)  combined;  03  <<  1,  I/O  <<  1. 

The  small  anisotropic  Rossby  number  limit  (Roa  <<  1)  covers  uniformly  all  four  cases 
(a)-(d).  The  limit  Roa  <<  1  is  not  a  mathematical  ethereal  limit  never  reached  in  “real 
physical”  geophysical  flows;  on  the  contrary,  our  rigorous  error  estimates  are  of  order  Roa 
where  the  anisotropic  Rossby  number  Roa  =  03/fl  =  a^Ro  is  of  order  10“^  —  10““*  for 
important  classes  of  atmospheric  and  oceanic  flows;  for  example,  03  w  10“^  and  Ro  w  0.1 
for  large  scale  atmospheric  flows  (synoptic  scales)  at  midlatitudes  which  implies  Roa  ~  10“^ 
(e.g.  Holton  [1992]). 

For  small  anisotropic  Rossby  number  flows,  we  rigorously  show  that  solutions  of  3D 
Euler/Navier-Stokes  equations  can  be  decomposed  as 

U(t,ii,a;2,a;3)  =  U(t,ii,i2)  +  V(t,xi,  12,3:3)  +  r 

where  U  is  a  solution  of  the  2D  Euler  system  with  vertically  averaged  initial  data;  r  is  a 
remainder,  its  Sobolev  norm  is  estimated  from  above  by  a  quantity  proportional  to  Roa,  for 
smooth  enough  initial  data. 

The  norm  of  V(t)  is  conserved,  ||V(t)||  =  ||V(0)|1  =  |1U(0)  —  U(0)|1,  therefore  it  is  not 
small.  The  vector  field  V(t)  is  related  to  v(i)  by  a  Van  der  Pol  type  transformation  based 
on  the  Poincare  propagator  and  v(t)  is  a  solution  of  a  new  system  which  we  caU  Extended 


Euler/Navier-Stokes  system.  The  Poincare  propagator  (Poincare  1910)  is  the  unitary  group 
solution  E(flt)  (E(0)  =  Id  is  the  identity)  to  the  Poincare  problem: 

di^  +  2fiJ#  =  - Vp,  V  •  $  =  0  (3) 

$(t)  =  E(-fit)$(0). 

We  use  Fourier  series  expansions  for  velocity  fields  U(x)  =  {U^{x),U'^{x),U^{x)) 

U(a:)  =  ^exp(t(niii  +  n2X2/a2  +  n3X3/a3))U„  (4) 

n 

where  are  the  Fourier  coefficients,  (ni, 722,713)  6  are  wavenumbers.  We  assume  that 
functions  have  zero  average  over  the  periodic  parallepiped. 

Poincare  [1910]  solutions  can  be  written  in  Fourier  space  as:‘ 

Eni^t)  =  {E(flt)}„  =  cos(-^2nt)I  +  sin(-^2fiORn  (5) 

a^\n\  |n|  0,3]'^] 

where  the  matrix  is  the  Fourier  transform  of  the  curl  operator;  we  denote  \fi\^  = 
n\  +  +  Tz^/a^.  The  mathematical  theory  of  the  Poincare  problem  (3)  has  attracted  a 

considerable  amount  of  attention  starting  from  the  work  of  Sobolev  [1954]. 

The  Poincare  velocity  v(^)  is  related  to  V(t)  by  V(^)  =  E(-~f2t)v(t),  which  is  a  canonical 
unitary  transformation  preserving  energy,  helicity  and  divergence  free  property.  For  fields 
which  depend  only  on  xi,  X2  but  not  on  X3,  E(fli)  =  Id.  Before  any  averaging  the  equations 
for  V  are  Euler-like,  but  with  time-dependent  coefficients.  The  extended  Euler  equations 
(which  are  obtained  after  time- averaging  for  non-resonant  domains)  are  for  the  velocity  fields 
which  are  three-dimensional  and  three-component  (3D-3C:  three  components,  dependence 
on  three  variables  xi,  X2,  X3).  The  extended  system  possesses  infinitely  many  conservation 
laws.  These  conserved  quantities  are  adiabatic  invariants  for  the  classical  3D  Euler  equations 
in  a  rotating  frame  in  the  small  anisotropic  Rossby  number  case.  If  the  vertically  averaged 
initial  data  are  reflectionally  symmetric  (that  is  satisfy  zero  flux  boundary  conditions  in  the 
half-sized  box),  we  can  fully  integrate  the  extended  system  and  obtain  the  exact  solutions 
for  v(t)  and  V(t).  We  underline  that  the  splitting  is  global,  the  initial  data  are  not  assumed 
to  be  small  and  strong  vertical  shearing  is  allowed. 
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As  usual  we  denote  X3-averaging  of  a  function  U  by  U,  and  for  vector  fields  U(xi ,  X2)  = 

prove  the  theorem  on  behavior  of  solutions  of  3D  Euler 
equations  which  gives  a  rigorous  sense  to  the  Taylor-Proudman  theorem  for  time-dependent 
flows.  Namely,  we  have 

Theorem  1.1. 

Let  02  be  non-resonant.  Let  U(t)  be  an  exact  solution  of  3D  Euler  equations  (1)  with  initial 
data  U(0)  =  U(0,xi,  X2,X3)  =  Uo(xi,X2,  X3).  Let  cr  -  9  >  a  >  3/2,  let  Mq^  >  0.  Then 
there  exists  Ti  >  0  such  that  if  ||U(0)||o.  <  Mq^,  then 

l|U(0-U(0ll"  <C’i2oa  ,  (6) 

for  0  <  ^  <  Ti;  C,Ti  depend  only  on  Moa-  Here  1|  •  \\o-  is  the  norm  in  a  Sobolev  space 
(of  divergence  free  periodic  functions)  and  U(i)  is  a  solution  of  2D-3C  Euler  equations 
with  vertically  averaged  initial  data  U(0)  =  U(0,a:i,a:2)  =  Uo(a:i,X2);  the  2D-3C  Euler 
equations  are  for  velocity  fields  independent  of  2:3: 

6,  U  +  U  •  VU  V  •  U  rr  0  (7) 

where  U  •  V  is  the  classical  2D  Euler  advection  operator.  In  particular, 

(8) 

dxi  dx2 


This  theorem  shows  that  the  exact  xa-averaged  vector  field  U(^)  is  close  to  a  solu¬ 
tion  U(i)  of  2D-3C  Euler  equations  (7)-(8)  with  vertically  averaged  initial  data  U(0)  = 
Uo(xi,  X2).  We  use  the  terminology  of  W.C.  Reynolds  and  Kassinos  [1994].  In  their  termi¬ 
nology  2D-3C  fields  have  three  components  and  depend  on  two  variables  x\  and  X2;  3D-3C 
fields  have  three  components  and  depend  on  three  variables  Xi,  X2  and  X3.  A  direct  corol¬ 
lary  of  this  theorem  is  separate  approximate  conservation  of  the  energy  of  components  U(t) 
and  U(t)  -  U(t)  of  the  flow.  This  implies,  in  particular,  that  U(t)  —  U(t)  is  not  small.  The 
existence  of  clatssical  conserved  quantities  for  2D  Euler  equations  (integrals  of  functions 
of  the  vertical  component  of  curl)  and  the  above  theorem  imply  existence  of  approximate 


conservation  laws  (adiabatic  invariants,  Arnold  1978)  for  3D  Euler  equations  in  the  case  of 
small  Roa- 

The  difference  U(t)  —  U(t)  (or  U(t)  —  U(<))  can  be  rigorously  approximated  by  V(t)  = 

E(— flt)v(t)  where  v(t)  is  a  solution  of  the  extended  Euler  equations.  These  equations  are 
for  the  velocity  fields  which  are  three-dimensional  and  three-component  (3D-3C).  They 
have  the  following  form,  where  U  is  the  solution  of  the  2D-3C  Euler  equations  (7)-(8)  with 
X3-  averaged  initial  data: 

atV  =  Be,(U,v);  (9) 

here  Bex(U,v)  is  a  bilinear  operator  and  v  =  0.  After  U  is  found  from  Eqs,  (7)-(8),  v 
is  obtained  from  Eq.  (9)  which  is  a  linear  equation  for  v  with  time  and  space  dependent 
variable  coefficients  determined  by  U. 

The  bilinear  operator  Bex(U,  v)  is  a  non-local  bilinear  operator  for  the  Z3-dependent  v 
field  and  U.  It  comes  from  a  non-local  restriction  of  Euler-like  operators  to  a  special  class 
of  wavenumber  interactions  determined  by  the  resonance  conditions.  These  Euler-like  oper¬ 
ators  were  explicitly  written  and  its  properties  discussed  in  Babin,  Mahalov  and  Nicolaenko 
[1995a].  Similar  extended  equations  were  written  for  the  Navier-Stokes  equations. 

The  remarkable  property  of  the  extended  system  (9)  is  that  (for  generic,  non-resonant 
a2,  03)  it  splits  into  an  infinite  sequence  of  independent  (uncoupled)  subsystems  of  linear  8- 
component  ordinary  differential  equations  (ODE^s)  for  the  corresponding  (divergence-free) 

Fourier  modes.  For  a  given  index  (711,712,713),  713  7^  0  they  describe  interactions  of  the  ^ 

quadruplet  of  Fourier  modes.  Each  subsystem  couples  8  non-autonomous 

ode’s.  The  equations  are  (in  the  case  of  non-resonant  domains) 

5iVn  =  i  [-Pn(Ufc  •  Th)Vm  +  X  ((A:  X  Ufc)  X  (m  X  v„i))],  (10) 

r7i|  =  t  ^2  ~  ^*^2 

m3  =  n3,  fc  -f  m  =  n 

where  n  =  (721,712/02,  7^3/03)  and  similarly  for  fc,  7h;  Pn  is  the  matrix  of  Leray  projection 
onto  the  space  of  divergence  free  periodic  vector  fields.  The  coefficients  of  this  system  are 
determined  by  the  solution  XJ(tjXi,X2)  of  2D-3C  Euler  equations  (7)-(8).  The  equations 
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are  very  simple  and  can  be  used  for  numerical  computations  of  3D  rotating  flows,  since  the 
subsystems  are  decoupled  in  a  very  convenient  way  for  parallel  computing. 

Every  subsystem  (10)  preserves  energy  and  helicity  separately,  which  gives  infinitely 
many  conservation  laws  for  the  whole  system;  moreover,  the  3D  .H^^-Sobolev  norms,  for 
every  s,  are  conserved  for  v  (including  enstrophy).  The  existence  of  infinite  number  of 
conservation  laws  for  the  extended  Euler  equations  is  in  contrast  to  the  classical  3D  Euler 
equations,  where  only  energy  and  helicity  are  conserved  (Serre  [1984a]  and  Serre  [1984b]). 

These  conserved  quantities  of  the  extended  Euler  equations  are  the  approximate  adiabatic 
invariants  (Arnold  1978)  of  3D  Euler  equations  in  the  small  anisotropic  Rossby  number 
situation. 

Being  so  elegant  and  even  explicitly  integrable  as  we  describe  below,  the  extended  equa¬ 
tions  describe  exact  3D  Euler  flows  with  high  accuracy  if  Roa  is  small.  We  show  that  U(f) 
can  be  rigorously  approximated  by  U(<)  +  V(t)  where  v(f)  is  a  solution  of  the  extended 
equations  (9),  (10).  The  vector  field  v(t)  is  related  to  V(t)  by  means  of  the  explicitly  given 
linear  unitary  operator  E(D<),  Eq.  (3)  (the  Poincare  propagator),  V(/.)  =  E(-flt)v(t). 

We  prove  the  following  theorem  delineating  the  structure  of  exact  solutions  of  3D  Euler 
equations  in  the  small  anisotropic  Rossby  number  situation;  effectively,  we  show  that  U  -1- 
V  U  as  Roa  — >•  0,  uniformly  with  respect  to  initial  data  and  strongly  in  Sobolev  norms, 
and  we  give  a  power  estimate  of  the  remainder: 

Theorem  1.2. 

■& 

Let  02)  ®3  be  non-resonant,  let  cr-o  >  35,  a  >  3/2,Mo<r  >  0.  Let  ||U(0)||<y  <  Mq^.  Let  U(t)  ''' 

be  an  exact  solution  of  3D  Euler  equations  (1)  with  U(0)  =  Uo(xi,X2)2:3)-  Let  E(fi<)  be 
the  Poincare  propagator,  v(t)  be  the  solution  of  the  extended  Euler  system  with  initial  data 
v(0)  =  U(0)  —  U(0).  Then  the  difference  U(t)  —  U(t)  —  E(— nt)v(t)  satisfies  the  estimate 

||U(0  -  U(t)  -  E(-nt)v(t)||2  <  CRoa  for  0  <  t  <  Tr  (11) 

where  Ti,  C  depend  on  Mq^- 


Similar  estimate  holds  for  solutions  of  3D  Navier-Stokes  and  extended  Navier-Stokes 
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equations  (uniformly  in  the  viscosity  u,  0  <  i/  <  1). 

Using  Theorem  1.2  on  approximation  of  solutions  of  3D  rotating  Euler  equations  by 
solutions  of  extended  Euler  equations  we  have  proven  regularity  of  solutions  of  3D  Euler 
equations  with  arbitrary  large  initial  data  on  arbitrary  long  time  intervals  in  the  small 
anisotropic  Rossby  number  situation: 

Theorem  1.3. 

Let  the  domain  be  nonresonant;  M  >  0,  T*  >  0  be  arbitrary  large.  Then  there  exists 
Rol  =  Rol{M,T*)  such  that  if  ||U(0)||38  <  M  and  0  <  oa/H  <  Ro^  there  exists  a  unique 
regular  solution  U(<),  0  <  <  <  T*  of  3D  Euler  equations  which  belongs  to  Hjs  as  0  <  t  <  T*. 
For  M  fixed,  T*  -»•  +oo  with  Ro*  ->  0.  Simultaneously  we  can  take  arbitrary  large  (but 
bounded)  sets  of  initial  data:  M  — ^  +oo  if  Ro"  —>■  0. 


The  problem  of  global  regularity  of  solutions  of  3D  Navier-Stokes  equations  has  been 
extensively  studied  by  many  mathematicians  and  still  is  an  outstanding  unsolved  problem 
of  the  applied  analysis.  Using  the  detailed  description  of  dynamics  of  the  extended  Navier- 
Stokes  equations  and  error  estimates  (analogue  of  Theorem  1.2  for  the  Navier-Stokes  equa¬ 
tions)  we  have  solved  the  problem  in  the  situation  of  a  small  anisotropic  Rossby  number  . 
Namely,  we  consider  the  3D  Navier-Stokes  equations  with  the  time-independent  forcing  F 
on  the  infinite  time  interval  : 

dtV  +  tycuTlcurlV  +  V-VV  +  2a3V  =  -  Vp -f- F,  (12) 

V-U  =  0. 


We  prove 
Theorem  1.4. 

Let  the  domain  be  nonresonant.  Let  F  G  M  >  0.  There  exists  Ro^{M^iyj\\F\\zs) 
such  that  if  llU(0)lli  <  M  and  0  <  aa/fi  <  then  there  exists  a  unique  regular  so¬ 
lution  U(0,  0  <  ^  <  +0O  of  3D  Navier-Stokes  equations  which  belongs  to  Hi  t  >  0, 
liU(t)lli  <  Ci(M,i/,11F1|38),  and  to  i/38  as  t  >  0;  and  1|U(0I|38  <  C'38(M,i/,  llFl|38,<o)  for 
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every  t  >  to  >  0.  (The  constants  Ci,  C^s  s^re  uniform  in  time  t,) 


In  Theorem  1.4  we  can  take  arbitrary  large  M,  ||F||38  and  a  solution  is  regular  if  Roa  = 
a^/n  is  small  enough;  but  the  question  whether  or  not  for  a  small  fixed  Roa  solutions  with 
arbitrary  large  initial  data  blow  up  in  Hi  in  finite  time  is  still  open.  In  all  cases,  the  energy 
estimate  ^  /J \\U{t)\\idt  <  t/^^\\V{0)\\‘^/T  +  holds. 

If  we  take  T  >  T"  =  liU(0)||g/i/,  =  2u-^\\F\\l^  +  2  then  |lU(t,)|li  <  Mi  for  some 

€  [0,r),  Then  applying  Theorem  1.4  with  a  shift  in  time  and  using  the  equality  of  a 
weak  solution  of  3D  Navier-Stokes  equations  with  the  regular  one  we  conclude  that 
a  weak  solution  with  arbitrary  large  initial  data  in  Hq  is  regular  fort  >  T*  and  is  attracted  to 
the  global  attractor  if  0  <  Roa  <  /io^(Mi,  r/,  HFHas).  The  regularization  time  T*  depends 
only  on  the  T2-norm  of  initial  data  and  the  critical  Rossby  number  Ro^  does  not  depend 
on  initial  data  at  all. 

We  do  not  need  to  impose  conditions  of  the  thin  domain  type  and  do  not  assume  that 
initial  data  are  close  to  2D  initial  data  in  some  Sobolev  space.  Furthermore,  the  dependence 
of  initial  velocity  fields  on  X3  can  be  arbitrary  and,  in  particular,  we  allow  strongly  sheared 
flows  in  2:3  (provided  the  fields  are  smooth).  The  parameter  a3  need  not  be  small;  flows 
corresponding  to  initial  conditions  may.  have  strong  dependence  (strong  X3-  shearing)  and 
large  -component. 

Remarkably,  if  the  vertically  averaged  initial  data  U(0)  are  reflectionally  symmetric, 
we  can  fully  integrate  the  extended  Euler  system.  Now  we  give  explicit  formulas  which 
express  solutions  of  the  extended  Euler  equations  (9),  (10)  in  terms  of  its  time-dependent 
coefficients.  The  corresponding  U(^,xi,X2)  +  V(i,  xj,  X2)  2:3)  are  exact  explicit  solutions  of 
3D  Euler  equations  in  the  small  anisotropic  Rossby  number  limit. 

Consider  the  case  when  X3-averaged  initial  data  have  the  reflection  symmetry  U^{—Xi)  = 
U^(xi)^i  ^  j,  f/*(— Xj)  =  =  1,2  at  t  =  0;  it  corresponds  to  the  usual  Euler 

zero-flux  boundary  condition  =  0  as  Xj  =  0  and  Xj  =  Taj.j  =  1,2  (which  are  the 
natural  boundary  conditions  for  the  2D  Euler  equations  in  a  rectangular  domain)  at  the 
boundary  of  the  half-sized  box.  Moreover,  we  take  U  (0)  odd  with  respect  to  the  reflection 
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symmetries  in  xi  and  12;  this  corresponds  to  the  no-slip  condition  for  the  third  averaged 
component.  For  (711,712,123)  we  denote  L{ni,n2,n3)  =  (-721,-722,723),  n'  =  (721,722,0). 

In  the  ca.se  of  symmetrical  averaged  data  the  extended  Euler  equations  decouple  into 
systems  of  ODE’s  describing  interactions  between  the  Fourier  coefficients  v„  and  VLn-  Let 


^n{t)  = -icurl\J)2^,  •  63 

be  the  Fourier  coefficient  of  the  vertical  component  of  curl  corresponding  to  the  wavenumber 
272'  =  (2721 , 2722, 0)  and 

^  _  n\  + filial  -  nl/al-3  .. 


Vnit)  = 


-uut) 


where  t^2n'(0  third  component  of  the  solution  of  2D-3C  Euler  system,  |n|^  =  nj  + 

n\la2  +  The  functions  ^  and  7/  are  real- valued. 

Then  Eqs.  (10)  split  into  equations  of  the  form 


(^n(0Rn  -  Vn{t)  ijnll') 


where  v„  =  =  (-fL, ^L)  and 


\  /  n  15  \  I  °  -n:ila3  722/03 

U  1  ,  U 

j  —  723/03  0  — 72l 

10/  \  ^  ®  / 

'  '  \  -722/02  72i  0 


Here  I  is  the  identity  matrix  and  zR.,1  is  the  curl  in  Fourier  space;  il'  and  are  skew- 
symmetric  matrices  which  commute.  The  commutativity  property  allows  us  to  write  general 
solutions  of  the  system  (13)  explicitly.  The  solution  is  given  by 


where 


-  ^2(0  *I«|I']) 


n(0  =  /  ^n{s)ds,  T2(t)  =  I  T]n{s)ds. 

Jq  Jo 


Vn(0) 

Vn(0) 
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We  note  that 


where  Id 


n{  +  n\la\  +  n^/aj;  we  always  take  divergence-free  v„.  For 


Navier-Stokes  equations  without  force,  the  right  hand  side  of  (15)  has  to  be  multiplied  by 
exp(-z/|npt). 

We  recall  that  V(t)  =  E(-flt)v(t).  An  important  observation  is  that  the  Poincare 
propagator  E(fit)  commutes  with  exp(^^^ri(t)R^)  and  exp(^^^r2(<)i|n|I')  and,  there¬ 
fore,  the  phases  2flt,  ri(t)  and  r2(<)  are  added  in  the  final  formulas  for  V(<).  Thus,  the 
vector  field  Y{t)  which  approximate  U(t)  -  U(t)  is  phase-locked  to  the  phases  2flt,  Ti{t) 
and  T2{t).  The  phases  T]{t)  and  r2(<)  are  associated  with  vertically  averaged  vertical  vortic- 
ity  curl\J{t)  ■  and  velocity  U  (t);  the  latter  is  multiplied  (in  Fourier  space)  by  the  wave 
operator  nf  -f  n|/a^  -  n§/a§  and  l/|n|  (smoothing).  The  phase  formulas  (17)  are  strikingly 
similar  to  the  Poincare  propagator  formula  (5);  except  that  they  describe  mode  coupling 
and  2Q,t  is  replaced  by  ri(t),  respectively  r2(<),  which  are  phases  associated  to  the  passive 
scalars  curlV{t)  ■  63,  respectively  U^{t).  The  latter  are  passively  advected  scalars  by  2D 
turbulence. 

Our  exact  solutions  completely  solve  the  outstanding  problem  of  parametrization  of  3-D 
rotation  dominated  turbulence  in  terms  of  pure  2-D  turbulence  of  verticaUy  averaged  fields, 
in  the  limit  of  small  anisotropic  Rossby  number  Roa-  Three-dimensional  rotating  turbulence 
in  the  small  anisotropic  Rossby  number  situation  decouples  into  irreducible  phase-locked 
turbulence  for  'V{t,xi,X2,X2)  and  2D  turbulence  for  vertically  averaged  fields.  Statistics 
of  2D-2C  turbulence  are  the  sole  player  for  U(t).  This  vindicates  all  current  research  on 
intermittencies,  vortex  dynamics  and  inverse  cascades  in  2-D  turbulence.  It  makes  very 
relevant  all  the  investigations  of  2D  vortex  dynamics  and  statistics  of  2D  tu'^bulence  by 
Kraichnan,  Montgomery,  SafFman,  Someria,  Zabusky,  McWilliams  among  many  others. 
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Experimental  studies,  direct  numerical  simulations  (DNS),  large-eddy  simulations  (LES), 
and  closure  approximations  have  established  that  solid  body  rotation  suppresses  the  non¬ 
linear  energy  cascade  from  large  scales  to  small  scales.  The  effect  of  rotation  is  through 
“phase  scrambling”  for  the  wave  phase.  For  MHD  turbulence,  Kraichnan  [1965]  pointed 
out  that  the  propagation  of  the  Alfvenic  fluctuations  disrupts  phase  relation  and -thereby 
may  be  expected,  on  the  average,  to  decrease  energy  transfer.  Similarly,  uniform  rotation 
causes  plane  waves  to  propagate  with  phase  speed  2fl&3/|^:|  resulting  in  “phase  scrambling” 
thanks  to  strong  dispersion  and,  on  the  average,  decreases  energy  transfer.  Our  results  on 
regularity  of  3D  Euler  and  Navier-Stokes  equations  reveal  the  regularizing  effect  of  rotation; 
the  mechanism  of  the  regularization  is  the  weakening  of  nonlinear  interactions  between  al¬ 
most  all  (non-resonant)  Fourier  modes  through  time-averaging.  The  mechanism  principally 
differs  from  previously  known  regularizing  mechanisms  where  viscosity  plays  the  crucial  role 
in  regularization.  The  theorems  proved  in  our  work  give  a  rigorous  mathematical  sense  to 
nonlinear  “phase  scrambling”. 

1.2  Inertial  Stability  of  Numerical  Algorithms  for  Fluid  Dynamics 

Exponential  attractors  give  a  better  answer  for  the  dynamical  and  numerical  significance 
of  long-time  flow  simulation  by  supercomputers.  Contrary  to  appearance,  neither  classical 
attractors,  nor  approximate  inertial  manifolds  are  uniformly  continuous  under  perturba¬ 
tions;  they  are  only  upper  semi-continuous.  As  a  consequence,  the  classical  attractor  of  a 
discretized  numerical  algorithm  or  nonlinear  Galerkin  scheme  cannot  uniformly  converge  to 
ti  a  exact  (continuum)  attractor  (in  the  Navier-Stokes  context).  In  contrast,  exponential 
attractors  are  robust  under  perturbations  (uniformly  continuous)  as  proven  for  Galerkin  ap¬ 
proximations  of  2D  Navier-Stokes  equation  in  Eden,  Foias,  Nicolaenko  and  Temam  (1994). 

The  existence  and  uniform  convergence  of  approximate  exponential  attractors  for  semi- 
implicit  and  explicit  discretization  schemes  for  the  Navier-Stokes  equations  are  of  particular 
interest. 

Exponential  attractors  are  enlargements  of  global  attractors;  they  are  in  general  strongly 
continuous  with  respect  to  perturbations,  wherea.s  the  latter  are  usually  only  upper-semicontinuous. 
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An  inertially  stable  algorithm  is  one  which  (i)  admits  a  discretized  exponential  attractor  and 
(ii)  the  latter  uniformly  converges  to  an  exact  exponential  attractor  for  infinite  resolution. 
An  inertially  stable  algorithm  shall  preserve  the  statistics  of  turbulence.  Inertial  stability 
of  an  algorithm  for  classical  2D  Navier-Stokes  equations  implies  the  same  for  smaU  Rossby 
number  rotating  3D. Navier-Stokes  equations.  Algorithms  implicit  for  the  Stokes  operator, 
yet  fully  explicit  for  the  Euler  advection  term  are  in  fact  globally  dissipative;  even  though 
energy  and  enstrophy  are  not  preserved.  Their  algorithmic  inertial  stability  condition  is 
superior  to  that  of  algorithms  semi-implicit  in  the  advection  term;  even  though  the  latter 
do  preserve  energy  and  are  trivially  globally  dissipative.  This  is  due  to  better  control  of 
turbulent  fluctuating  direct  cascades  of  enstrophy  and  inverse  cascades  of  energy  by  the  first 
class  of  algorithms.  Optimal  inertial  stability  yields  effective  control  of  unresolved  subgrid 
scales,  and  optimal  resolution  of  large  scale  vortices  in  3D  small  Rossby  number  rotating 
turbulence. 

VVe  have  investigated  “inertial  stability”  of  numerical  algorithms  as  the  proper  tool 
to  ensure  the  uniform  approximation  of  dynamical  turbulent  attractors.  Our  substantial 
results  (Babin  and  Nicolaenko  [1994],  Babin,  Mahalov  and  Nicolaenko  [1996b])  within  the 
context  of  two-dimensional  Navier-Stokes  turbulence  shall  provide  a  handsome  pay-off  in 
the  context  of  geophysical  turbulence,  through  the  exact  2-D,  2-C  “geostrophic”  operators. 
We  shall  expand  this  effort. 

Let  us  contrast  “inertial  algorithmic”  stability  versus  classical  stability  concepts  for 
numerical  schemes. 

For  numerical  accuracy,  one  traditionally  considers  the  error  estimates  for  trajectories. 
Exponential  growth  in  time  is  generally  present  in  the  clctssical  results.  As  shown  in  the 
most  recent  and  sharpest  estimates  for  the  errors  induced  by  both  classical  and  nonlinear 
Galerkin  methods  on  N.S.  flows,  (Devulder-Marion-Titi  [1994],  Rautmann  [1980]),  a  non¬ 
linear  Galerkin  trajectory  will  shadow  an  exact  one  for  N.S.  flows  up  to  a  finite  time,  for 
a  given  tolerance.  Of  course,  the  sharper  the  algorithm,  the  longer  the  shadowing  time. 
Such  a  separation  is  to  be  expected:  even  when  there  are  no  computational  errors  at  all, 
generally  a  small  perturbation  of  the  initial  condition  can  cause  exponential  splitting  of  two 
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trajectories. 

Much  work  has  been  done  on  errors  in  estimating  the  global  dynamical  attractor  /I, 
for  an  extensive  variety  of  approximation  schemes;  they  are  all  based  on  the  upper  semi¬ 
continuity  of  the  global  attractor. 

Exponential  attractors  give  a  better  answer  for  the  dynamical  and  numerical  significance 
of  long-time  flow  simulation  by  supercomputers.  They  are  continuous  with  respect  to  the 
full  Hausdorff  (Federer,  [1969])  set-distance,  for  “good”  approximation  schemes. 

The  existence  of  exponential  attractors  for  turbulent  Navier-Stokes  flows  has  a  major  impact 
on  their  large  scale  computing.  For  instance,  with  spectral  methods  it  is  common  that  half 
of  the  modes  (the  small  ones)  carry  less  than  to  10”^  of  the  total  energy.  Hence 

the  computer  tends  to  be  saturated  by  small  wavelengths  carrying  little  energy.  Similarly, 
with  finite  differences,  the  unknown  varies  little  between  two  successive  grid  points  when 
Dx  ~  10“^.  Exponential  attractors  filter,  in  large-scale  computing,  these  subgrid  scale 
effects  built  into  the  discrete  system,  A  numerical  algorithm  must  not  only  be  vested  with 
conventional  stability;  as  a  discretized  large  finite-dimensional  dynamical  system,  it  must 
have  built-in  uniform  stability  of  its  own  discretized  exponential  attractor,  VVe  implement 
such  a  program  based  on  the  entirely  new 

Definition.  An  approximate  exponential  attractor  is  the  exponential  attractor  of  an  ‘Hner- 
tially  stable^  scheme. 

Definition.  A  scheme  is  Snertially  stable^  if  and  only  if: 

(i)  it  is  stable,  dissipative  and  possesses  an  exponential  attractor  Mi^  of  its  own.  Here 
N  is  the  dimension  of  the  phase  space  of  the  scheme. 

(ii)  for  all  e,  there  exists  N{e)  so  that: 

dH{Mi^,M)  <  e 

{Hausdorff  distance). 

iii)  the  constant  of  exponential  convergence  rate  to  {for  dynamics  of  the  approximate 
scheme)  is  of  the  same  order  as  estimates  of  that  constant  for  the  exact  theoretical  M. 

In  ii),  €  is  the  tolerance,  iii)  ensures  that  approximate  dynamics  lock  onto  the  approxi¬ 
mate  M/v  os  fast  as  the  exact  dynamics  onto  M. 


Corollary.  After  time  transients  of  the  same  order^  the  approximate  and  exact  dynamics 
take  place  on  Mj^  and  M. 

Remark.  There  is  no  need  to  further  compute  Myv-  After  uniform  transients^  the  dynamics 
of  the  scheme  are  on  Mat.  This  is  not  true  for  convergence  to  A  and/or  Aj^. 

Eden,  Foias,  Nicolaenko  and  Temam  [1994]  have  proven  that  the  classical  Galerkin 
scheme  is  “inertially  stable”  (this  is  a  semi-discrete  scheme  where  the  time  variable  is  kept 
continuous).  This  is  the  litmus  test  for  any  other  scheme: 

•  Determine  which  approximation  schemes  are  vested  with  an  approximate 
exponential  attractor  and  are  indeed  inertially  stable? 

The  task  is  much  less  formidable  than  it  appears,  since  the  difficult  questions  of  finite- 
time  stability  and  dissipativity  of  whole  classes  of  algorithms  have  already  been  resolved  in 
the  literature  quoted  above. 

Our  approach  is  a  novel  one  for  the  effective  control  of  unresolved  subgrid  scales  in  terms 
of  long  time  integration  of  turbulence.  As  a  payoff,  statistics  of  turbulence  are  uniformly 
preserved  for  long-time  integration.  Also,  at  any  time  (especially  long  times  in  the  DNS), 
there  is  a  finite  piece  of  an  exact  physical  trajectory  tracked  by  the  computed  one.  Among 
concrete  and  more  practical  new  results  (Babin  and  Nicolaenko  [1994],  Babin,  Mahalov, 
Nicolaenko  [1996b])  we  find  that: 

•  inertial  stability  of  a  given  algorithm  crucially  depends  on  the  built-in  spatio-temporal 
physical  scales  (including  unresolved  subgrid  scales)  of  turbulence  for  the  specific  physical 
model  computed  (usually  not  true  for  classical  algorithmic  stability).  For  classical  (incom¬ 
pressible)  Navier-Stokes,  one  finds  that  the  physical  scales  of  the  inertial  range  of  turbulence, 
as  well  as  the  extent  of  that  range,  play  a  key  role.  Numerical  turbulent  attractor  stability 
must  be  dealt  with  on  a  basis  of  specific  classes  of  models:  the  physical  scales  of  a  turbulent 
boundary  layer  around  an  aircraft  wing  span  differ  from  those  of  quasi  2-dimensional  (2-D) 
geophysical  flows. 

One  also  finds  that  (at  least  for  incompressible  Navier-Stokes): 

•  the  better  the  numerical  scheme  approximates  all  conservation  laws  of  the  underlying 
inviscid  Euler  equation,  the  more  performing  it  is  from  the  turbulent  attractor  stability 
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perspective.  One  conjectures  this  to  be  true  for  more  general  algorithmic  models  and  plan 
to  test  it  for  geophysical  models.  For  instance,  semi-implicit  (in  the  nonlinear  advective 
term)  schemes  for  Navier-Stokes  do  conserve  energy  exactly,  but  fail  to  preserve  ‘‘enstrophy” 
(that  is,  the  total  energy  of  the  vorticity  field).  They  turn  out  to  be  much  worse  off  than 
fully  explicit  (in  the  nonlinear  advection  term)  schemes  which  approximately  yet  uniformly 
preserve  both  energy  and  enstrophy. 
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